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Abstract  -  Echoes  reflected  from  the  sea  floor  hold 
a  great  deal  of  information  beyond  what  is  visible  in 
simple  time  series  or  spectra.  In  addition  to  reflections 
from  different  boundary  layers,  structures  contained  in 
these  layers  will  be  excited  by  energy  from  the  incident 
pulse  and  may  reradiate  this  energy.  This  reradiated 
energy  will  contain  information  about  the  physical  na¬ 
ture  of  the  structure  and  surrounding  media.  This 
work  looks  at  application  of  time-scale  signal  expan¬ 
sions  to  extracting  such  information  from  fathometer 
echoes. 


1  Model  for  Air  Bubbles  in 
Bottom  Sediments 

A  commonly  occuring  sea  bottom  structure  is  consti¬ 
tuted  by  gas  bubbles  suspended  in  sediment.  An  in¬ 
cident  sound  wave  on  a  gas  bubble  surrounded  by  a 
fluid  medium  will  excite  the  bubble;  that  is,  induce 
oscillations  in  the  bubble  at  the  constituent  frequen¬ 
cies  of  the  sound  wave.  The  degree  of  excitation  of  the 
bubble  is  dependent  on  the  relationship  between  the 
spectral  content  of  the  incident  wave  and  the  natu¬ 
ral  frequencies  of  the  system  defined  by  the  immersed 
bubble.  Energy  imparted  to  the  bubble  will  be  maxi¬ 
mum  at  frequencies  corresponding  to  those  natural,  or 
resonance  frequencies. 

If  the  product  of  the  wavenumber,  ko,  at  the  bub¬ 
ble’s  fundamental  resonance  and  bubble  reudius,  r,  is 
such  that  kor  <<  1,  the  immersed  bubble  oscillating 
at  its  fundamental  resonance  frequency  can  be  mod¬ 
elled  as  a  simple  one  degree  of  freedom  mass-spring- 
damper  system  ([Devin  59],  [Med win  76],  (Junger  86] 
pp.  67-69)  with  impulse  response  described  by  the 
linear  second  order  differential  equation 

Af  C  +  aujMC  +  ETC  =  ^(0  (1) 

E^^uation  1  describes  the  motion  of  the  system  in  re¬ 
sponse  to  an  impulse  excitation  6{t),  at  frequency  u;, 

^This  work  was  supported  by  NRL  under  the  MTEDS 
program. 


with  M  the  system  mass,  K  the  stiffness  of  the  system, 
and  aujM  describing  the  energy  dissipation  within  the 
system.  The  term  a  is  called  the  system  damping 
constant.  This  model  assumes  the  energy  dissipation 
is  proportional  to  the  bubble  volume  velocity,  The 
solution  to  Eq.  1  is  a  simple  exponentially  damped 
sinusoid 

C(t)  =  (2) 

with  tt(t)  the  step  function,  /  =  u;/27r,  and  ujd  = 
vA  —  oP-jAw.  Note  for  a  <<  I,  Ud  ^  u,  implying 
that  for  light  damping  the  resonance  frequency  of  a 
submerged  bubble  with  radius  r  can  be  approximated 

r  _  V^jPoP  ,,, 

Here  p  is  defined  as  the  additive  density  of  the  ele¬ 
ments  comprising  the  medium  around  the  bubble,  7 
is  the  ratio  of  specific  heats  for  the  bubble  gas,  and 
Po  is  the  static  pressure  in  the  medium.  Following  the 
development  in  [Devin  59],  the  effective  system  mass 
is  due  to  the  medium  surrounding  the  bubble  surface 
(this  assumes  that  the  mass  of  the  gas  inside  the  bub¬ 
ble  is  negligible)  and  is  given  by 


(4) 


The  stiffness  associated  with  the  system  is  due  to  the 
compressibility  of  the  gas  that  forms  the  bubble  and 
is  described  by 


„  _  3-irPo 
4irr3 

The  system  damping  constant,  a,  relates  the  energy 
dissipation  mechanisms  at  work  on  the  system.  For 
volume  pulsations  of  submerged  bubbles,  there  are 
three  dominant  mechanisms  for  energy  dissipation. 
These  are  thermal  damping  due  to  heat  conduction 
between  the  bubble  gas  and  surrounding  liquid,  ra¬ 
diation  damping  as  energy  is  transported  away  from 
the  bubble  acoustically,  and  viscous  damping  result¬ 
ing  from  viscous  forces  at  the  bubble-fluid  junction. 
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Detailed  discussions  of  the  damping  mechanisms  and 
expressions  for  a  can  be  found  in  [Devin  59]  and 
[Medwin  76).  System  damping  constant  will  be  a  sum 
of  damping  constants  due  to  each  of  the  three  mecha¬ 
nisms  cited  above 


Q  =  Qt  +  Q!r  +  ctv  (6) 

where  Qt,  Or,  and  Ou  refer  to  damping  due  to  thermal, 
radiation,  and  viscous  mechanisms,  respectively. 

The  significance  of  the  bubble  impulse  response  de¬ 
scribed  in  Eq.  2  is  that  the  center  frequency  and  decay 
rate  are  characteristic  of  the  geometry  of  the  bubble 
and  the  physical  properties  of  the  bubble  and  medium. 
Such  information  can  be  exploited  for  sediment  anal¬ 
ysis  applications. 


2  Wavelet  Techniques 

The  continuous  wavelet  transform  (CWT)  is  a  tech¬ 
nique  which  expands  a  signal  onto  a  time-scale  plane 
via  correlation  with  a  set  of  basis  functions,  called 
wavelets,  indexed  by  time,  6,  and  scale,  a.  Simply  put, 
a  wavelet  can  be  any  function,  /i(i),  which  is  bandpass 
and  zero  mean.  Families  of  wavelets  are  formed  by 
dilating,  or  scaJing-,  and  shifting  h{t)  in  time 

Coefficients  are  generated  at  each  point  a,  6  in  scale 
and  time  by  correlating  with  a  wavelet  that  has  been 
scaled  and  shifted.  For  a  signal  s{t)  the  CWT  will 
consist  of  a  set  of  coefficients,  (t>a,h{s)y  generated  by 
taking  an  inner  product  of  the  signal  and  wavelet  at 
every  point  in  scale  and  time  with  (l>a,b(s)  given  by 

CSO 

</>a,6(s)  =  /  s{t}hl^b{t)dt  (8) 


In  continuous  form,  the  CWT  is  an  energy  preserving 
transform  and  will  yield  true  time-scale  energy  distri¬ 
butions  for  a  signal.  Discrete  implementations  of  the 
CWT  will  yield  good  approximations  of  a  signals  time- 
scale  energy  content  if  the  time-scale  plane  is  sampled 
densely  enough  [Daubechies  90).  Note  also  that  if  the 
spectral  content  of  the  wavelet  used  is  unimodal  with 
energy  that  is  well  localized  in  frequency  there  is  an 
inverse  relationship  between  scale  and  frequency.  This 
property  allows  CWT’s  generated  by  such  wavelets 
to  be  treated  as  time-frequency  expansions.  Because 
of  their  ability  to  simultaneously  decompose  a  signal 
into  temporal  and  spectral  components,  wavelet  ex¬ 
pansions  atre  effective  in  non-stationary  signal  analy¬ 
sis. 

The  optimal  (with  respect  to  signal-to-noise  ratio} 
set  of  basis  functions  on  which  to  expand  a  signal  are 
ones  which  precisely  match  the  signal.  However,  in 
many  cases  this  is  not  possible  as  the  structure  of  a 


signal  cannot  be  completely  determined.  Nevertheless, 
basis  functions  that  closely  mimic  the  time-scale  be¬ 
havior  of  a  process  can  be  very  effective  at  localizing 
signal  component  energies.  Wavelet  transforms  offer 
the  ability  to  perform  isometric  signal  expansions  with 
a  flexibility  in  the  choice  of  the  basis  function.  This 
flexibility  can  be  exploited  to  yield  processing  gains 
and  unique  feature  extraction  capabilities. 

Processes  which  generate  transient  signals  in  the 
form  of  exponentially  damped  sinusoids  can  be  de¬ 
scribed  as  time-scale  processes  if  they  exhibit  damping 
constants  that  are  linearly  proportional  to  frequency. 
This  label  implies  that  signals  arising  from  these  pro¬ 
cesses  will  consist  of  components,  or  modes,  which 
possess  constant  time-bandwidth  products.  However, 
it  is  not  necessary  that  a  process  be  time-scale  in 
order  to  be  exploitable  by  wavelet  transform  tech¬ 
niques.  The  most  straightforward  method  of  tailoring 
a  wavelet  to  match  the  signal’s  time-scale  behavior 
is  to  use  a  wavelet  which  is  a  copy  of  the  signal.  In 
practical  applications,  a  quadrature  sampling  scheme 
will  be  necessary  to  deal  with  phase  differences  be¬ 
tween  the  sampled  signals  and  the  discrete  wavelets. 
For  processes  whose  components  generate  signals  in 
the  form  of  damped  sinusoids,  this  suggests  use  of  a 
damped  complex  exponential  with  added  terms  to  give 
a  wavelet  of  the  form 

g-(a/2+i)u,ty(jj  ^ 

Operating  at  ranges  where  q  1  and  a;  1  will 
make  the  contributions  of  the  last  two  terms  in  eq.  9 
negligible  and  allow  the  wavelet  to  be  approximated 
as 

-  6)  (10) 

y/a 

Now  note  that  for  light  damping,  the  impulse  re¬ 
sponse  of  submerged  gas  bubbles  given  in  Eq.  2  can 
be  rewritten 

C(t)  =  •  (11) 

The  scaling  operation  of  the  wavelet  expansion  de¬ 
scribed  in  the  previous  section  produces  basis  func¬ 
tions  in  which  damping  is  linearly  proportional  to 
frequency.  This  implies  that  wavelet  expansions  us¬ 
ing  damped  complex  exponential  wavelets  are  optimal 
with  respect  to  signal-to-noise  ratio  for  processes  pro¬ 
ducing  damped  sinusoidal  component  signals  that  ex¬ 
hibit  a  constant  damping/frequency  ratio;  i.e.  pro¬ 
cesses  that  are  time-scale  in  nature.  In  the  insteince 
of  bubbles  suspended  in  fluid  sediments,  this  requires 
that  damping  at  resonance  is  linearly  proportional  to 
the  resonance  frequency.  For  light  damping  this  re¬ 
quirement  means  that  damping  at  resonance  is  in¬ 
versely  proportional  to  bubble  radius,  i.e.  that  a  is 
constant.  Under  these  circumstances,  bubble  radius 
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wiU  be  proportional  to  the  wavelet  scale  variable.  Sys¬ 
tem  damping  for  submerged  gas  bubbles  is  not  lin¬ 
early  proportional  to  resonance  frequency;  however, 
over  certain  bandwidths,it  may  be  reasonably  approx¬ 
imated  as  linearly  proportional  to  resonance  frequency. 

Figure  1  is  the  time  series  and  power  spectrum  of 
a  fathometer  echo.  Visible  in  the  time  series  is  the 
penetration  of  multiple  layers  on  the  sea  floor.  The 
spectrum  indicates  a  broadband  return.  In  order  to 
examine  features  representative  of  the  sea  bottom  it 
will  be  necessary  to  deconvolve  the  fathometer  source 
signal  from  the  return.  Such  deconvolved  echoes  will 
characterize  the  impulse  response,  or  transfer  function 
of  the  sea  bottom.  Shown  in  Fig.  2  is  the  wavelet 
transform  of  the  return  shown  in  Fig.  1  deconvolved 
from  the  source  signal.  This  transform  used  the  ba¬ 
sis  functions  of  Eq.  10  to  match  the  bubble  impulse 
responses  described  in  Eq.  11.  Note  the  numerous 
clusters  of  high  energy  coefficients,  or  modes,  in  the 
signal.  These  modes  are  indicative  of  concentrations 
of  resonators  buried  throughout  the  bottom.  Here,  the 
resonators  are  gas  bubbles  submerged  in  the  bottom 
sediment  which  have  been  excited  by  the  incident  fath¬ 
ometer  pulse.  Deconvolving  the  source  signal  yields 
the  bubble  impulse  responses. 

Wavelet  transforms  yield  many  important  clues  use¬ 
ful  for  identification  and  analysis  of  bottom  sediments. 
Resonators  such  as  submerged  gas  bubbles  can  be 
identified  and  localized  in  time  (depth)  and  frequency. 
Because  of  the  basis  function  used  and  the  nature 
of  the  bubble  impulse  response,  modal  decay  rates 
can  be  estimated  from  wavelet  transform  coefficients. 
Modal  decay  rates  can  be  exploited  as  discriminators 
of  the  presence  of  submerged  gas  bubbles.  Addition- 
adly,  anaJysis  of  center  frequencies  and  decay  rates  may 
be  useful  in  estimating  bubble  sizes.  Wavelet  trans¬ 
form  coefficients  can  also  be  used  to  estimate  the  en¬ 
ergy  in  each  mode  of  a  signal. 

For  damped  complex  exponential  wavelets  applied 
to  damped  sinusoidal  signals,  examination  of  the  tem¬ 
poral  behavior  of  wavelet  transform  coefficients  allows 
estimation  of  system  damping,  a.  Assume  a  damped 
sinusoidal  signal  s{t)  which  starts  at  time  with  am¬ 
plitude  A,  damping  constant  a,,  and  center  frequency 
It  can  be  demonstrated  that  after  peaking,  the 
wavelet  transform  coefficients  will  decay  at  a  rate  that 
is  solely  a  function  of  the  signal  decay  rate, 

[Priebe  95].  For  a  particular  mode,  the  decay  rate 
can  be  determined  by  observing  the  temporal  vector 
of  wavelet  transform  coefficients  containing  the  modal 
peak  as  it  decays  in  magnitude  beyond  that  peak  (Fig. 
3).  If  <j>a,hi  and  0a,  62  are  two  such  coefficients  from 
that  vector  separated  in  time  by  6t  —  a/{b2  —  bi)  then 
Ois  can  be  estimated  by 


_  2a  |0a,62l  a 

'  |0a.6j(62-6i) 


(12) 


Gas  bubbles  trapped  in  bottom  sediments  often  occur 
in  large  concentrations  that  are  excited  concurrently, 


giving  rise  to  coexisting  multiple  modes.  Consider 
individual  signal  mode  with  center  frequency 
rate  a^,  and  start  time  which  is  surrounded^n  ti 
and  frequency  by  other  modes  with  amplitudes  A 
center  frequencies  decay  rates  a„,  and  start  tirn^ 
tn .  The  energy  from  these  nearby  modes  will  influ 
ence  the  wavelet  transform  coefficients  correspondin 
to  the  mode  at  u;,,  The  degree  of  this  influence  wiU 
be  proportional  to  the  modal  separations  in  time  and  ^ 
frequency.  When  evaluating  modal  decay  rates  for  sig-  - 
nals  with  closely  spaced  modes  (such  as  the  signal  of 
Fig.  2),  the  estimates  obtained  from  the  wavelet  trans¬ 
form  coefficients  should  be  treated  as  moving  averages 
in  time  and  scale. 

Figure  4  plots  the  modal  damping  coefficients  (q) 
as  a  function  of  frequency  for  the  modes  of  the  sig¬ 
nal  of  Fig.  2,  a  return  from  sediment  containing  sub¬ 
merged  gas  bubbles.  Analysis  of  returns  of  known  bub¬ 
bly  sediments  gives  a  good  indication  of  what  range  of 
modal  decay  rates  to  expect  for  these  systems.  Also 
on  the  plot  of  Fig.  4  is  a  curve  showing  the  theoret¬ 
ical  values  for  damping  constant,  q  as  a  function  of 
resonance  frequency.  The  theoretical  curve  was  gen¬ 
erated  bas^xl  on  expressions  developed  in  [Devin  59]. 

It  assumed  constant  depth  air  bubbles  submerged  in 
sediment  whose  density  was  twice  that  of  water.  Fur¬ 
ther,  the  curve  also  assumed  spherical  bubbles  oscil¬ 
lating  at  their  fundamental  resonance  frequency.  Note 
that  there  is  about  an  order  of  magnitude  discrep¬ 
ancy  between  the  theoretical  curve  and  the  mean  of 
the  measured  values.  This  discrepancy  may  be  due 
several  factors,  some  of  which  are:  the  bubbles  were 
generally  not  spherical  in  shape  [Abegg  94]  ,  they  were 
distributed  over  a  range  of  depths  (radiation  damping 
increases  significantly  with  depth),  the  gas  involved 
was  typically  methane  (the  curve  assumed  air),  and 
the  sediment  density  differed  from  the  value  used  for 
the  curve. 

Shown  in  Fig.  5  are  two  probability  density  func¬ 
tions  (pdf’s)  for  modal  decay  rates.  They  are  based 
on  measured  modal  decay  rates  from  two  known  bot¬ 
tom  types  “  hard  surface  and  gassy  sediment.  Signal 
modes  connected  with  reflection  at  interfaces  will  ex¬ 
hibit  higher  decay  rates  which 'tend  to  be  more  char¬ 
acteristic  of  the  source  transfer  function  (as  evidenced 
by  the  pdf  corresponding  to  the  hard  surface  data). 
Decay  rates  from  the  gassy  sediments  tended  to  be 
more  lightly  damped.  No*te  that  the  pdf’s  did  not  dis¬ 
tinguish  modal  sources  within  the  individual  echoes. 
All  modes  present  within  a  ping  were  used  to  generate 
the  pdf  for  that  particular  bottom  type.  Hence,  more 
highly  damped  modes  resulting  from  reflection  at  in¬ 
terfaces  are  included  in  the  pdf  for  gassy  sediment. 

Figure  5  was  used  to  implement  a  simple  bottom 
classification  scheme.  By  analyzing  fathometer  echoes 
from  a  ship  moving  at  about  4  knots  and  sounding  once 
a  second  a  bottom  map  of  a  section  of  Eckernfoerde 
Bay  was  generated.  This  map  is  displayed  in  Fig.  6. 
The  map  was  formed  by  averaging  modal  decay  rates 
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from  deconvolved  fathometer  echoes  over  celk  of  0.1 
meter  resolution  in  depth  within  each  echo.  Any  cells 
with  three  or  more  modes  present  whose  average  decay 
rate  fell  between  0.005  and  0.0328  were  classified  as 
gassy  sediments  and  are  shown  as  gray  areas  on  the 
map.  Cells  with  three  or  more  modes  present  whose 
average  decay  rate  exceeded  0.0328  were  classified  as 
interfaces.  These  are  the  white  areas  on  the  map. 

Other  useful  features  can  be  extracted  from  fath¬ 
ometer  echo  wavelet  transforms.  These  include  modal 
density  (simply  the  number  of  modes  in  an  echo)  and 
penetration  (defined  as  the  maximum  range  difference 
between  modes  in  an  echo).  Pdf’s  for  each  of  these 
features  are  shown  in  Figures  7  and  8  for  echoes  from 
hard  surfaces  and  gassy  sediments.  While  modal  den¬ 
sity  and  penetration  may  not  be  unique  to  any  par¬ 
ticular  bottom  class,  using  these  two  features  along 
with  modal  decay  rates  may  provide  useful  sea  floor 
classification  tools. 


[Abegg  94] 


3  Conclusions 

The  wavelet  transform  offers  a  practical  tool  for  echo 
sounder  signal  analysis  because  of  its  ability  to  decom¬ 
pose  these  signals  into  their  constituent  modes  and 
extract  relevant  features  from  those  modes.  By  exam¬ 
ining  the  modal  features,  judgments  can  be  made  as 
to  the  nature  of  any  underlying  structure  in  the  sea 
floor.  These  methods  have  been  demonstrated  to  be 
particularly  suited  to  detecting  submerged  gas  bub¬ 
bles.  Further,  use  of  multiple  features  available  from 
echo-sounder  signal  decomposition  has  potential  for 
wider  applications  to  sea  floor  classification. 
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Figure  1 :  Fathometer  echo  from  gassy  sediment. 


Figure  2:  Wavelet  transform  of  fathometer  echo  of  Figure  L  Coefficient  energy 
levels  are  mapped  to  colors.  Signal  modes  are  clearly  visible. 


Figure  3:  Estimation  of  modal  decay  rate  of  a  damped  sinusoidal  signals  by  examination  of  the  temporal 
vector  of  wavelet  transform  coefficients  containing  the  modal  peak.  This  technique  assumes  use  of  a  damped 
complex  exponential  wavelet. 
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Figure  4:  Plot  of  damping  constants  for  signal  modes  of  impulse  response  shown  in  Figure  2  vs.  bubble 
resonance  frequencies.  Overlaid  on  the  plot  is  a  curve  showing  theoretical  damping  constants  for  air  bubbles 
sediment  at  20m  depth  for  band  of  interest.  Below  the  frequency  scale  is  a  scale  showing  the  bubble  radius 
95sociated  with  the  resonant  frequencies  of  the  theoretical  curve. 
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Figure  5:  Probability  density  functions  for  modal 
decay  rates  evident  in  echoes  from  known 
hard  surfaces  and  gassy  sediments. 


range  (m) 

Figure  6;  Seafloor  classification  map  of  modal  decay  rates  for 
section  of  Eckernfoerde  Bay. 
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Figure  7:  Probability  density  functions  for  modal  densities  evident  in  fathometer  echoes  from  known  hard 
surfaces  and  gassy  sediments. 


Figure  8:  Probability  density  functions  for  penetration  evident  in  fathometer  echoes  from  known  hard  surfaces 
and  gassy  sediments. 
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